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Abstract 



The effective properties of piezoelectric laminates have been analyzed, based on 
the calculation of internal fields and making use of a simple matrix manipulation 
method. The results are expressed in a compact notation which is convenient for 
numerical implementation and at the same time suitable for further analytical treat- 
ments. A detailed analysis of fully compatible ferroelectric domain structures shows, 
that the results for arbitrary piezoelectric laminates can be further simplified and 
specific property relationships for rank-1 laminates of tetragonal and rhombohedral 
crystals are derived. The method is finally applied to the analysis of various hierar- 
chical domain structures. Detailed orientation relationships between the particular 
domains in some important domain pattern are given to make these structures 
accessible for the presented method. Some numerical results for tetragonal barium 
titanate illustrate the effects of different domain arrangements on the effective prop- 
erties. 
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1 Introduction 



Modern functional materials, e.g. ferroelectrics, are often characterized by an 
internal microstructure, which in many cases can be described as a laminar 
structure. Although this is usually a simplification which neglects some effects, 
it is a useful approach that allows an analytical analysis of the overall proper- 
ties of these materials. It provides insight into the physical properties of the 
composite, which may strongly differ from the local, single-phase properties. 



In particular, we are interested in the linear material properties of a laminar 
compound and consider elastic, dielectric and piezoelectric coefficients as well 
as spontaneous strain and polarization. The general approach may be applied 
also to other coupled problems and associated linear properties, if the relevant 
boundary conditions are appropriately chosen. In this context, magneto elec- 
tric laminar composites (jAvellaneda and Harshel . Il994 iNan et all 1200.4 e.g.) 
should be mentioned 



The relation between microstructure, local and effective prop erties is one of 
the central subjects of the theory of composites (jMilton . Laminates or 

stratified materials are the simplest conceivable composites with a variation 
of the material properties in only one direction, the direction of lamination 



n fjMiltonl . 120021) . These simple type of com posites are also called rank-one 



laminates (iMiltonl . l2002h or 2-2 composites ( Newnham et all Il978t ). For an 
overview of the historical development of the lamination theory and general so - 
lution strategies for various coupled problems see, for instance. iMilton (|2002h . 



The knowledge of the effective behavior of piezoelectric multilayers is a crucial 
problem not only for stratified piezoelectric materials, which have been in 



vestigated for various t r ansducer applicat i ons dCao et al. L . 9931: Zhang et al. 
1994 I.Tanas and Safari Il99.4 iLous et all 1200(1 iKoh et all 1200,4 . A soecia 



special 

case of a laminar piezoelectric structure is a ferroelectric crystal with a lami- 
nar domain structure, i.e. a ferroelectric solid, which consists of two different 
ferroelectric domain types (same crystal structure in the two phases, but with 
different crystal orientations with respect to the interface) separated by a 
crystallographically defined interface. The domains are arranged in laminar 
configuration and separated by plane interfaces (domain walls). An example 
of such a structure is shown in fig. 1. The picture also demonstrates, that real 
microstructures are often much more complex. However, laminar stuctures 
can be used as building blocks in the description of more realistic patterns 
by assembling differently oriented laminates into a complex, often hierarchical 
microstructure of a single grain or crystallite. 



In a recent review, Bhattacharva and Ravichandran ( 2003| ) have pointed out 
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the importance of the knowledge of the overall behavior of ferroelectric do- 
main patterns for large strain piezoelectrics. The problem of laminar ferro- 
electric domain structures becomes also an important research object in the 



context of engineer e d domain structures dPark et al. , 199S ; Wada et al. L 199S ; 



Zhang et all 1200 ll lYin and Caol. bnni iLiu and LiL booJ 12004 lid and Liu 



2004 iDavis et allbnnd IWada et all 12005b . 




Fig. 1. SEM image of a ferroelectric ceramic (tetragonal PZT), showing a domain 
structure, which consists of differently oriented blocks of laminar domain patterns. 

The purpose of this contribution is to provide a self-contained presentation 
of the theory of piezoelectric laminates, ranging from arbitrary stratified and 
electromechanically coupled multi-layers to simple as well as hierarchical fer- 
roelectric domain patterns. The theory offers both effective material proper- 
ties (elastic, piezoelectric and dielectric constants) and average spontaneous 
strains and polarizations as well as internal fields within the layers and allows 
an analysis of effective (or averaged) symmetry properties of a multi-layer or 
multi-domain crystal. We do not restrict ourself to certain crystal symmetries 
or volume fractions. 



Vario us var i ants of the problem and solutions have been published i n the 



past (E urik L 197dt Turik and Bondarenkol. 19741 Aksakava and Farnelll . 1988 ; 
Kreherl. 119981 lOibianskv and Torouatol Il999t lErhart and Caol Il999l 12001 : 
Bednarcvkl . 120031 ILiu and Lil . 120031 : lid and Liul . 120041 ) . Although the solution 



tion 

of internal fields in piezoelectric multi-layers and the calculation of effective 
prope rties is straigh tforward and can be analytically derived without any diffi- 
culty ( Milton . 2002| ) . this problem is worth to be revised again. In particular, 
because the analytical expressions are generally quite intricate, any contri- 
bution which provides a simplified approach should be welcome. Additionally, 
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the various authors concentrate on various but separate aspects of the problem 
and not all published results are correct. 



The author has used the presented mod e ls in previous works on tetragonal 
ferroelectric ceramics ( Rodel and Kreher , 2003| ). Here, we generalize the ap- 
proach to cover a variety of cases, ranging from arbitrary piezoelectric multi- 
layers to particular ferroelectric domain structures with different crystal sym- 
metries. 



After stating the problem in section 2, it is therefore useful to review in section 
3 a solution for the general problem of a piezoelectric multilayer with arbitrary 
number and structure of the constituent layers, which may apply for arbitrary 
piezoelectrically coupled laminates. The aim of our approach is to notate the 
internal fields and effective properties in a compact form, which is, at the same 
time, easy to implement numerically and convenient for further analytical 
treatments. 



In section 4 we will consider the degenerate case of a fully compatible fer- 
roelectric domain structure where certain orientation relationships are met. 
Further, we will give some examples for special cases with particular crystal 
structures. Then, domain patterns with incompatible polarizations are ana- 
lyzed in contrast to the fully compatible domain structures. Finally, we will 
briefly consider hierarchical laminates, in particular some important cases of 
hierarchical domain structures. 



In the case of ferroelectric domain structures, one has to note, that domain 
walls are generally mobile and can move under applied fields. This may cause 
additional non-linear contributions to the macroscopic behavior, which is be- 
yond the scope of this paper. However, the proposed approach provides also 
a convenient basis to describe non-linear reversible and irreversible response 
due to reversible or irreversible domain wall motion. 
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2 Problem formulation and notation 



We consider a piezoelectric solid (see fig. 2) under so-called homogeneous (and 
quasistatic) boundary conditions, i.e. boundary condition at the outer bound- 
ary of the solid are chosen so that internal fields of stress, cr(x), strain 7(2;), 
electric field E(x) and dielectric displacement D(x) would be homogeneous, 
if the material properties would be homogeneous. 



t X 3 




Fig. 2. Piezoelectric body subjected to prescribed tractions and electric potentials. 
The body has a laminar micro-structure. 

The electric fields and the linear elastic strain are defined by the gradients of 
the electric potential <j)(x) and elastic displacements u(x), respectively: 



_d_ 

dx-. 
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Surface tractions, t, may be applied at the outer boundary with normal vector 
n. The tractions are related to the elastic stress stress tensor <r(x) 

ti{x) = «7y(aj)nj(x) (2) 



Stress fields and electric fields have to fulfill the equilibrium conditions 

=0, ^O(-) = 0, (3) 
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where we neglect volume forces and free charges. The equilibrium conditions 
imply the following equation at internal interfaces: 



ni[taj|] = 0, ni[|A|] = 0, (4) 

where [\X\] denotes the jump of the quantity X at the interface and n is the 
normal vector of the interface plane. 

We define the macroscopic fields by the volume averages over the local fields: 



a V 



- J tr(x) dx = (tr(x)) , 7 = 7/ / 7 ^ dx = W x )) 

V V 

E = ^j E(x) dx = (E(x)) , D =^J D ( x ) dx = W x )) > ( 5 ) 



V V 

which are the so-called micro-macro-relations. 



The local fields are related to each other by the linear piezoelectric constitutive 
law: 



7y = £yki °"ki + d i]k E k + 7;j (6. a) 

D i = d iU a kl + e° k E k + P* (6.b) 

where is the elastic stiffness tensor under constant electric fields, e? k is 
the dielectric permittivity tensor under constant stress and d-^ k is the tensor 
of the piezoelectric strain constants, relating the strain to the electric field. 



Ferroelectric crystals exhibit a spontaneous strain and spontaneous polariza- 
tion, 7 s und P s respectively, at temperatures below the curie temperature. 
The spontaneous strain is sometimes called stress free strain or eigenstrain, 
because it is caused by the ferroelectric phase transformation only, and not 
by any applied stress. The reference state is a cubic state. Correspondingly, 
the spontaneous polarization is the polarization of a crystal without applied 
electric fields. In the case, that thermoelastic and pyroelectric effects are of 
interest, the spontaneous strain and the spontaneous polarization can be ex- 
tended using the thermal expansion and pyroelectric coefficients. 



One should be careful if other field quantities will be kept constant. For in- 
stance, the strain under constant dielectric displacement is not equal to the 
that one under constant electric field. 



The constitutive equations can be conveniently expressed using Voigt's short 
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index notation (reduced index notation). A compact notation of the electrome- 
chanical constitutive law, is obtained by combining equivalent mechanical and 
electrical field quantities in vectors with 9 components: 

















q s 




q = 













/ 7 s 
/a 



q = Q p + q s , Q 

and the inverted constitutive law is 



f CE J . 

°a,/3 u a,] 



p = R (q - q s 



R = Q" 1 



/ r D -h 



a = 1,2,. ..6 



i = 1,2,3 



a, (3 = 1,2,. ..6 
i, j = 1,2,3 



a,P= 1,2, ...6 
i,j = 1,2,3 



(7) 



(9) 



Combined property matrices, Q and R are symmetric. Here, we have com- 
bined strain and dielectric displacement, as well as stress and electric field, 
because they are thermodynamically equivalent and we will see later, that 
some results will be significantly simplified, if stress and electric fields are 
prescribed. However, it should be noted that mechanical stress and dielectric 
displacement as well as strain and negative electric field are equivalent from 
mathematical point of view: 





7a Id 
-Ef\ 



(10) 



Lf +j s 



e i/3 _£ ?j 



The above introduced combination of field quantities is applied to local fields 
as well as to the macroscopic, averaged fields. The intention of this work is to 
calculate the effective material properties which fulfill the following equation 
describing the macroscopic linear material behavior: 



q = Q* p + q* 



;i2) 



with 



P 



(P(*)> 




q 




<q(*)> 



(7(*)> 
\{D(x)) 
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q* s = , Q 




Following an often used ansatz in the homogenization theory we will express 
the local fields by disturbance of the external field 

p = Ap + p A , (13) 

where A is an interaction matrix, which describes the interaction of the various 
phases under applied loads p and p A is a vector which describes stresses and 
electric fields, which arise in the lamella (k) due to incompatible spontaneous 
stresses and incompatible spontaneous polarizations. 

The general formalism to calculate the effective material properties is to in- 
clude (8) and (13) into (12) 



q = (q) = (QP + q*> = (Q( Ap + p A ) + q s ) = (QA) p + (Qp + q s ) , (14) 

from which we derive the general expression for the effective properties by 
comparison with equation (12): 

Q* = (QA>, q* s = (Q P A + q s ) (15) 



The remaining task is to calculate the coefficients of the matrix A and the 
vector p A 
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3 Arbitrary piezoelectric laminar structures 

In order to derive a solution for a piezoelectric multilayer further assump- 
tions are necessary. We consider a piezoelectric solid, which has a layered mi- 
crostructure consisting of n individual layers. Each layer k may have different 
electromechanical properties: 
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Inside each layer, the material properties are constant. 

The interfaces between the layers are considered to be flat, plane and parallel 
to each other. As the direction of lamination we have chosen the ^-direction, 
i.e. the macroscopical ^-direction should be normal to the interfaces. Addi- 
tionally, we assume that every layer is very thin compared to its length. Then, 
we can neglect effects at the edges and the volume fraction of of the kth layer, 
q k , is proportional to it's thickness. 

Under this assumptions, we can presume that the internal fields do not vary 
in the directions perpendicular to direction of lamination, i.e. the internal 
fields depend only on x 3 . Then, the equilibrium equations imply that the 
following field components are constant within the layers and continuous at 
the boundary between two layers i and j, which means, that they are equal to 
the applied fields: 
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Since the electric field has to be curl free and strains have to be compatible, the 
following fields components are continuous across the interfaces and following 
conditions must be satisfied: 



y u = 7ii — qr = qr (i6.e) 

7 22 = 7 22 — q { 2 ] = q { 2 ] (16-f) 

7i 2 =7{ 2 — ??> = ? 6 W (16-g) 

= E\ = E 1 — ► p [ 7 k) = p 7 (16.h) 

E\ = E\ = E 2 — ^ fe) =p 8 (16.i) 



Additionally, within a particular layer, all fields are constant. As a result, the 
volume average of a property X is the sum of X^ k \ weighted with the volume 
fraction £w 

(X) = 



We summarize the conditions, which have to be satisfied: 

qW = qP i = 1,2,6,9 (17.a) 
p! fc) =Pi 1 = 3,4,5,7,8 (17.b) 

As we want solve a set of equations (13) for the internal fields p, we need four 
additional equations. We make use of the micro- macro-relation (5): 

(p? ) )=ih 1 = 1,2,6,9 (17.c) 



It is furthermore convenient to split the indices into two groups, according to 
which boundary condition holds for pf 

i = 1,2,6,9 
i = 3,4,5,7,8 

The calculation of the components of A and p A is straightforward and given 
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in detail in the appendix. We summarize the solution in the following form: 
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where we made use of the following abbreviations: Q is a 4 x 4 submatrix of 
the property matrix Q, and R is the inverse of Q. 

/ n (p) n (p) n (p) n (p) \ 
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and Q is the inverse matrix of the volume average of R: 



- \ -1 



Q = ( R) or in index notation Q-^ = (R~ 



? (p) 



Note that Q-^ is a 4 x 4 sub-matrix of the 9x9 matrix Q^. R~^ is the inverse 
of this 4x4 matrix Q~^, and not a 4 x 4 sub-matrix of Q _1 . R is symmetric: 
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If we introduce the interaction matrices into (15) we find 



Qk = 


/ \ 










= 




% = 








«a = 









ik 



nk/ ' X^iTl^Ik 



9f 



ij j 



+ (Of, 



ik 



(19.a) 
(19.b) 
(19.c) 

(19.d) 
(19.e) 

(19.f) 



This equ a tion f or th e effective p r opert ies are equivalent to those given by 
Liu and Li ( 2003f ). and Li and Liul ( 2004 ). Although their solution is very el- 



egant, our derivation offers some advantages, as we have directly given the 
equations for the internal fields p^. Additionally, we have also derived the 
the relationship for the effective sponta neous strain and polari zation of the 
composite. In the often cited paper by Erhart and Cad ( 1999j ) the authors 
have used wrong conditions at the interfaces, which obviously led them to 
wrong conclusions. The rather detailed treatment of the general case may be 
also seen as a correction of their results. 



Alternative sets of effective properties, like L* may be calculated with very 
similar formulas. The only difference is, that the separation of indices accord- 
ing to the boundary conditions, or likewise, the separatio n into sub-matrices 



as shown in the appendix, will be different (jKrehe 



yparatior 
3,|l993). 



To close this section, we remind, that the effective properties have to obey the 
simple bounds (Voigt-Reuss-type bounds) 

(Qr 1 < (QT 1 = R* < (R) (20) 



or 



W 1 < Q* < (Q) (21) 

where the inequalities are ment in the sense of positive definiteness. 
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4 Simple compatible domain structures 



Ferroelectric domains are a result of symmetry breaking phase transforma- 
tions. In the ferroelectric phase exist several equivalent variants with the 
same, but rotated, crystal structure. Spontaneous polarization and sponta- 
neous strain due to the phase transformation are strictly coupled to each 
other and correlated to the crystal lattice. The term ferroelectric refers to the 
fact, that the polarizations, and thus the strain, can be switched from one 
direction to an other by applied electric fields. 



For sake of energy minimization, ferroelectric crystals usually exist in a multi- 
domain state and the interfaces between the domains are also correlated to the 
crystal lattice. This yields particular crystallographic orientation relationships 
between different domains. 



Consider a material with arbitrary symmetry. On a microscopic scale (within 
the domains) the material is assumed to obey linear piezoelectric behavior 
with the single crystal properties Q and q s . We do not consider nonlinear 
effects, neither switching nor any kind of domain wall motions. It should be 
emphasized that the present treatment can easily extended to take into ac- 
count reversible and irreversible motions of the domain walls. This will be 
published elsewhere. 



In this section we want to consider a fully compatible laminar arrangement of 
two domain types. The term compatible generally means, that the spontaneous 
strain and polarization of the two domains are compatible at the interface so 
that there appear neither stresses nor electric fields inside the domains, if there 
is no macroscopic stress and no macroscopic field. It is obvious, that this is 
the case, if the the components 7^, 75^ and 7^ of the spontaneous strain and 
P| of the polarization are continuos across the interface, i.e. these components 
have to be homogeneous. This requirement can also be demonstrated based 
on eqn (18.e), which says that p = if all q? are homogeneous for i = 1,2,6,9; 
i.e. q? 1 = q? n ;i = 1,2,6,9. 



Assume, that we have a crystal in arbitrary orientation with respect to to 
the interface, with n = (0, 0, 1)', on one side of the interface. For a particu- 
lar crystal symmetry, only a limited set of orientations relationships between 
the crystal orientation and the interface (or domain wall) are possible. For 
a detailed analysis of the this relationships, including a refined compatibil- 
ity condition which takes into account the change of crystal symmetry due 
to the ferroelectric phase tra nsformation, further arguments are necessary 
( Shu and Bhattacharval . l200lh . However, this is beyond the scope of this pa- 
per. We simply assume, that the plane, characterized by n = (0, 0, 1)*, is a 
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Fig. 3. Orientation relationship of compatible ferroelectric domains. Thick arrows 
represent polarizations. Two possible types of unit cells indicate spontaneous strains. 

possible domain interface for the particular crystal orientation. 



A fully compatible domain can be created by a rotation of the crystal on 
the opposite side of the interface. The part of the crystal, which is in the 
original orientation, will be called domain I, denoted by the suffix I. likewise 
the rotated part will be domain II. The crystal coordinate system of domain 
II is then rotated by the following relationship: 

x lI = u>x l , (22) 

where a; is a rotation matrix 



-1 











-1 











1 



(23) 



It should be noted, that this rotation ensures only compatibility of spontaneous 
strain and polarization. It does not necessarily describe a pair of ferroelectric 
domains, which form as a result of a phase transformation. This question 
should be considered separately for particular crystal symmetries. 



15 



The properties of the rotated crystal with respect to the global coordinate 
system are expressed by 

p. 11 = c^p. 1 , or P n = ujP 1 

7J 1 = u ik Wji 7k b or 7 11 = u 7 1 w* 

^ijk — ^il ^jm ^kn rfj mn 
"ijkl = ^im ^jn ^ko ^lp S mnop 



We find the following relationship for vectors, like polarizations: 
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For symmetric 2nd rank tensors like strain and permittivity holds 
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which shows, that the proposed general orientation relationship between the 
two domains ensures compatibility of spontaneous strain and polarization, 
independent of the particular symmetry of the crystal. For the piezoelectric 
coefficients and the elastic constant we get (in Voigt's notation): 
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If we combine the properties of the domain II in the compact form we get 
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In other words: the rotation uj does not change the magnitude of the compo- 
nents of Q and q s , but it changes the signs of some components: 
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i,j = 1..9 

i = 1,2,6,9; j = 4,5,7, 
i = 3; j = 4,5,7, 

i = 1..9 



i = 4,5,7, 



We note furthermore, that the submatrix Q-j with i,j = 1,2,6,9 is homoge- 
neous, i.e. according to the classification in the appendix, the domain system 
is partially homogenous. As it has been shown in A, internal fields and ef- 
fective properties are considerably simplified in this case. It follows from the 
condition of compatibility of the q s that 







(24) 



Then, (19.e)-(19.f) yield 

gf* = <?f>, i = 1-9 or q s * = (q s ) 



(25) 



It is convenient the to separate the interaction matrix A^ p ^ = I + A'^ with I 
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beeing the 9x9 identity matrix and 



) 



i = 1,2,6,9; 



k 



3,4,5,7,8 (26) 



4 p) = 

ik 

A'i p) = 



i,k = 1,2,6,9; 



i = 3,4,5, 7,8; k = 1..9 



(27) 



(28) 



ik 



(29) 



We may write for the internal fields 



p (p) = p + A'p 



In this representation the internal fields are expressed by a homogeneous con- 
tribution (corresponds to the generalized Reuss approximation) and an addi- 
tional contribution A'^ p ^ caused by the material inhomogeneity (anisotropy) 

Additionally, the following components Q~ are homogeneous: 



and this yields A~~ = A n = A 23 = A 63 = A g3 = 0, and 
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It should be noted, that the interaction matrix A for a ferroelectric domain 
structure is an expression of the anisotropy of the linear material properties 
and of the crystallographic symmetry relationship between the two domains. 

We can express the effective properties of a domain structure by the formula 




Q* = (Q) + A 



(30) 
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with 



A Ik = 



A~ = 

ik 



ij J mk 



1=1,2,6,9 (31) 
k = 1..9 

1 = 3,4,5,7,8; (32) 
k= 1,2,6,9 

i,k = 3,4,5,7,8 (33) 



As the Qis = with i=l,2,6,9, are homogeneous, also the components 
A3; = Aj3 = 0. The symmetric matrix A;k has the following form: 
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Therefore, in a ferroelectric domain structure with arbitrary symmetry, the 
only components of the effective linear properties, which differ from the simple 
average (Q) (lower bound) are Q44 = £44, Q45 = S45, Q47 = d u , Q i8 = d 2 4, 
Q55 — 5*55, Q57 = c?i5, Qsg = g?25, Q77 — Q78 = ^12 and Qss = £22- 

The constants Sn, S12, S13, Siq, S22, S23, S33, S26, 5*66, ^31? ^32, ^36 an d ^33? 
of the single crystal, are homogenous and the effective constants are equal to 
the single crystal values. 



The effective constants S* i: S* 5 , S^, S^, S^ 5 , S^, S^g, S^ e , cfj^, d^i, d\ 2 , 
j* j* j* j* j* j* j* j* 

u 22' u 13> u 23> u 33> u 34> u 35> u 16> u 26 

fraction and are zero at £ = 0.5. 



^22 > ^*3> ^23' ^33' ^34; ^35; ^*6' ^26' e i3 an d £ 23 depend linearly upon the volume 



It should be emphasized, that the simplifications for ferroelectric domain struc- 
tures, which have been derived in this section, are not possible for alternative 
sets of effective constants, like R* or L*. The internal field components, which 
have to be continuous across the interface are only for applied stresses and 
electric fields equal to the macroscopic fields. Alternative set of constants have 
to be calculated using formulas similar to those derived in sect. 3. 
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In the following we want to study important cases for particular crystal sym- 
metries. 



4-1 Tetragonal 90° domain structures 



Barium titanate (BaTiOs), lead titanate (PbTiOa) and titantium rich lead 
zirconate titanate (PbZr x Ti!_ x 03, PZT) are important tetragonal piezo- and 
ferroelectric material, which often show laminar 90° domain structures, as 
sketched in fig. 3. 



The properties of the single domain single crystal are usually described in a 
crystal coordinate system with the tetragonal axis (pseud o-cubic [001] = po- 
larization direction) parallel to macroscopic £3 l|NveLll964 - We will denote the 
single crystal properties described in the crystal coordinate system by the suf- 
fix °: Q°, q so A tetragonal ferroelectric phase is characterized by spontaneous 
polarization and strain 
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There are six possible tetragonal variants, characterized by the polarization 
directions: [001], [100], [010] and their negative directions, and three sponta- 
neous strain states associated with the tetragonal phase transformation. We 
may distinguish two principle domain types: 90° domains, where the angle 
between the polarization directions is approximately 90°, the strain states are 
different and the domain walls are the pseudo-cubic {110} planes. The second 
type of tetragonal domains are 180° domains, which have identical strain in 
both domains, but the polarizations point in opposite directions. All planes 
with a normal vector perpendicular to the polarization are possible domain 
walls. This example is considered in sect. 4.2. 



We assume, that the single crystal in the pseudo-cubic parelectric state is 
rotated by 45° about the X2 axis, so that the pseudo-cubic (101) becomes 
parallel to the plane with n = (0, 0, 1)* There are four possible arrangements 
of tetragonal domains with this plane as a domain wall (figure 4). 



For the purpose of this work, it is sufficient to consider only the first example. 
The second one is easily obtained by interchanging the volume fractions of 
domain I and II, and the remaining two arrangements are the same as the 
first two, but with inverse polarization directions (rotated by 180°). 
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Fig. 4. The four possible 90° domain structures with a (101) | |n = (0,0, 1)* domain 
wall. 



In order to end up with a compatible spontane ous strain, the domain I is 
rotated the angle 6 (jShu and Bhattacharval . l200lh 



7T 
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with 



arccos 



2 7i°i7l3 
,7n +7n 



As the spontaneous strain is small the rotation <p is small and has negligible ef- 
fects on the effective properties. The orientation of domain II may be obtained 
from domain I by applying the orientation relationship (22) and (23). 

The, we find 7^, 7g 6 = and P| = in both domains. For the rotated linear 
properties holds e l2 , £23, d u , d 16l d 2 \, d 22l d 23l d 25l d 34 , d 36 , S u , S 16 , £24, 5*26, 
S34, S3Q, S45, = in both domains. 



This results in further components of the interaction matrix A to be zero: 
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Then follows that for tetragonal 90° domains the interaction matrix A^ p ^ has 
the form: 
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4-2 Tetragonal 180° domain structures 



A particular simple case of an 180° domain structure is obtained, assuming 
that the pseudo-cubic (100) plane is the domain wall. However, we want to 
consider the general case, where all planes with n _L P are allowed (see fig. 5 
a). Then, an additional rotation ip about the [001]-axis of the crystal should 
be taken into account and the transformation matrix, u/, which rotates the 
reference orientation into an arbitrarily oriented domain I, is: 



x 



UJ X 



CO 
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sin ip cos ip 
- cos if) sin ip 



/ 



Then, we find the crystal properties of domain I with respect to the macro- 
scopic coordinate system in terms of the single crystal properties in crystal 
reference coordinates, S°, d°, e°, 7 s0 , and P so : 
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S — [2(S n S 12 ) S G6 )} 
S' 22 = 5*12 — (cos 4 ij> — cos 2 ip) S' 
S' 23 = 5°! + (cos 4 ip — cos 2 ip) S' 

S' 2A = ~sm(^) S' 

S44 = 5*66 — 4 (cos 4 tp — cos 2 ip) S' 



In the case of a 180° domain it is interesting to note, that applying the orien- 
tation relationship (22) and (23), will yield two anti-parallel oriented domains. 
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Fig. 5. 180° domains, a) Interfaces between compatible, anti-parallel domains are 
not necessarily parallel to simple crystallographic planes, b) Compatible bi-crystal 
with tilted crystal lattices but anti-parallel polarizations, c) 180° domain wall with 
anti-parallel crystal lattices vectors, but no tilt. 



However, the resulting structure (fig. 5 b) is generally not a 180° domain struc- 
ture in the common sense, because the crystal lattices of the two domains are 
tilted with respect to each other. Although this does not violate the condition 
of compatibility of strain and polarization, we should denote this example 
rather as a bi-crystal then a domain structure. A domain structure with anti- 
parallel polarizations, which forms as a result of a symmetry breaking phase 
transformation (fig. 5 c) can be described by following transformation matri- 
ces: 

/ n n , \ 
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u>"x c 
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-1 



sin ip cos ip 



cos ip sin ip 







It is obvious, that the elastic and dielectric properties of the two domains are 
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identical, but all piezoelectric coefficients will have opposite signs. R is again 
homogeneous. 



We find, that the following components of the interaction matrix A become 
zero: A u , A 24 , A 64 , A m , A 15 , A 25 , A 65 , A 67 , A 97 , A 18 , A 28 , A 98 = 0, and the re- 
maining components are Ag 5 , An, A27, A 68 7^ 0. 
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The following components of A are zero: A 44 , A 45 , A 47 , A 48 , A 57 , A 58 and A 78 . 
The matrix A takes the form: 
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As a result, the effective properties of a 180° domain structure is very close 
to the generalized Reuss approximation. The three contributions A 55 , A 77 and 
A 88 can easily be explained. The piezoelectric shear effect will cause opposed 
dielectric displacements, if a macroscopic shear stress <t 5 is applied. This results 
in a reduced compliance of the composite. Anti directed in-plane shear of the 
two domains in the case of applied E 2 due to the inverse piezoelectric effect 
results in in-plane shear stresses which restrict the free deformation of the 
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stack. This results in a decrease of effective e^- If the volume fractions are 
equal, the effective dielectric constant in 2-direction, , equals the single 
crystal dielectric constant under fixed strain ■ A similar effect influences 
the dielectric constant in polarization direction, s\ x , where opposed elongation 
and compression in the two domains are caused by an applied E\. This effect is, 
of course, also present for other domain orientations, but most easily explained 
in the 180° case. 

For an 180° bi-crystal (fig. 5 b), S24 and £34 will have opposite signs. Then, the 
components A14 and A24 will not be zero. In this case, an applied shear stress 
04 will cause inhomogeneous internal stresses o\ and 02. Moreover, A44 and A47 
will be non-zero because due to the piezoelectric coupling, the inhomogeneous 
stresses 0\ and a 2 will result in dielectric displacements D\. This is particular 
interesting, as it means, that there may appear an additional piezoelectric 
coefficient d\ A for the case of the bi-crystal, although the coefficients of the 
uncoupled constituents are zero. 



25 



4-3 Rhombohedral domains 



Within the framework described above, it is easy to analyze other crystal 
symmetries. In particular, the treatment of rhombohedral domain structures is 
very similar to the tetragonal case. Rhombohedral phases of space group R3m 
and R3c appear in many solid solutions of lead titanate, like PbZr x Tii_ x 03 
with x > 0.52, or PMN-xPT (Pb(Mg 1/3 Nb 2 /3)xTii_ x 3 ) with x > 0.7. The 
rhombohedral phase is characterized by a spontaneous polarization along the 
pseudo-cubic [111] direction and a spontaneous strain which results in the 
rhombohedral distortion of the unit cell: 
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There exist two types of non-180 domain patterns, the first one, has the (001) 
plane as the domain interface and the angle between the polarization vectors 
P sI and — P sI1 is arccos | ~ 70.53°. In a second type of domain structures 
the (110) plane is the domain interface and the angle between P sI and — P sI1 
is arccos f — |) ~ 109.47°. This principal patterns may are sketched in (110) 
plane in fig. 6. 




[110] 

Fig. 6. Simple non-180° domain patterns in rhombohedral phases 



Due to the rhombohedral distortion, a small additional rotation by angle 
about the pseudo-cubic axis [110] is necessary to ensure strain compatibility. 
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For the case of a (001) domain interface <fi is 

7n + 7i 3 



arccos 



7fi 2 + 2 7 fi7?3 + 37f 3 



s 2 



and for the (110) interface <fi is 
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77!i 2 + 27! 
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Again, as the spontaneous strain is small the rotation 
negligible effects on the effective properties. 



is small and has 



The single crystal properties of a rhombohedral material are usually described 
in the so-called hexagonal setting, where the pseudo-cubic [111] axis is the 
crystal c-axis and p aralle l to th e macroscopic x$ direction and the crystal has 
trigonal symmetry (jNvel ( 1964 ). fig. 7). With respect to this orientation, the 
linear property matrix of a rhombohedral single crystal, Q h and qr sh is: 
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First, we want to consider a 70.53° domain structure with (001) interface. In 
order to apply the formalism of sect 4 we have to rotate the unit cell by —90° 
about the X3 axis and by an angle arccos (1/V3) w 55° about the X2 axis, 
yielding the orientation of domain I with the pseudo-cubic [110] direction 
parallel to X\. 
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Fig. 7. Perovskite unit cell in hexagonal (trigonal) setting, seen along the 
pseudo-cubic [111] direction (hexagonal c axis). The corresponding hexagonal unit 
cell is indicated by the lattice vectors a h and b h . The orientation of the pseudo-cubic 
cell is indicated by pseudo-cubic lattice vectors [100], [010] and [001]. 

Similarily, in order to describe the 109° domain structure with the (110) in- 
terface the rotation oj" will rotate the domain I so that n = (0, 0, 1)* is the 
interface normal to the pseudo-cubic (110) plane. 
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(34) 



Alternatively, the same result is obtained by subsequent rotations of the do- 
main I in the 71° setting by -90° and 180° about the X2- and X3-axis, respec- 
tively. Then, orientation relationship (22) and (23) can be applied in order 
to get the orientation and properties of domain II for both cases. The effec- 
tive properties of the rhombohedral domain pattern can be calculated using 
formulas of sect. 4. 

Note, that the chosen orientation of domain I is in principle a monoclinic or 
othorhombic setting of the unit cell, i.e. the reference orientation of mono- 
clinic or othorhombic crystals can be directly used to calculate the effective 
properties of monoclinic domain structures. 
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5 Charged domain walls 



Although, compatible domain structures are considered to be favored for sake 
of energy minimization, the existence of not fully compatible structures are 
quite common. For instance by applying sufficiently high electric fields to a 
single crystal, the domains can arrange themselves so that the spontaneous 
strains are compatible but the polarizations are not. This is possible, if the 
average polarization of all domains, which are favored by the direction of the 
applied field, is aligned with the field, but there is no possibility to arrange 
them so that the polarizations are compatible. For instance, in a tetragonal 
crystal, an applied electric field in [lll]-direction would favor the [100], [010], 
and [001]-domain, but there is no possibility to arrange this polarizations 
in a way, that gives continuous polarizations across all domain boundaries. 
However, the strains can be compatible. The concept of domain engineering 
makes use of this principle. We will study this particular example in detail in 
sect. ??. In this section we want to emphasize the differences between a fully 
compatible rank-1 laminate (sect. 4) and a rank-1 laminate with incompatible 
polarizations for arbitrary material symmetry. 




Fig. 8. Orientation relationship of compatible ferroelectric domains 



Consider a domain I, arbitrarily oriented with respect a plane with normal 
n = (0,0, 1)*, as in sect. 4 (fig. 8). A domain II with compatible strains can 
be created by rotating the material on the other side of the interface by 180° 
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about the xi-axis. Then, the rotation oj 
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(35) 



yield a mechanically compatible, but electrically incompatible domain struc- 
ture, if the strain component 7f 2 = 0. This means, in order to describe com- 
patible spontaneous strains, we have to chose the direction of the coordinate 
system so that the in-plane shear vanishes, which should be always possible. 



We obtain the following properties of domain II: 
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The domain wall is charged by discontinuous polarization components normal 
to the interface. 

p = n[\P s \] = 2Pi, 

which gives rise to internal electric fields which will ensure the continuity of 
dielectric displacements, even if the external field is zero. The internal fields in 
both domains are directed in opposite direction to ensure a macroscopic field 
free state. The internal electric field may also result in piezoelectric strains 
and generate dielectric displacements parallel to the domain wall. 



Obviously, the sub-matrices Qik, i, k — 1, 2, 6, 9 are inhomogenous for general 
piezoelectric materials. Then, the simplifications for partially homogeneous 
laminates, made in sect. 4 do not apply and the general solution from sect. 3 
has to be used in order to compute effective properties and internal fields. 
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6 Hierarchical domain patterns (rank-two laminates) 

Ferroelectric domains are often arranged in a more complex way. However, as 
compatibility of strains and dielectric displacements and energy minimization 
leads to essentially plane domain interfaces, the laminar structures play an 
prominent role. More complex patterns can be build from laminar structures 
by combining differently oriented laminates to so called hierarchical structures 
or rank-two- microstructures. This structures can be considered as laminates 
of laminates, i.e. laminates, in which each lamellae again consists of a fine 
laminar structure. 

Here we want to consider only microstructures which have compatible spon- 
taneous strains everywhere. We will denote the larger scale lamellae as bands. 
The interfaces between the bands are again flat. The numerical calculation 
of effective properties can be performed without difficulty following the lines 
drawn in the preceding sections. The remaining task is to define the orienta- 
tions of the individual domains and domain bands in an appropriate way. One 
should be aware, that the domain bands show a reduced symmetry compared 
to the individual domains. Therefore, we should expect even more reduced 
symmetries in the linear properties of the hierarchical composite. 

We want to consider some important examples of such types of structures. The 
first example is a combination of tetragonal 180° and 90° domains. We assume, 
that the fine structure is a 180° domain stack and two bands of 180° domains 
are arranged in such a way, that they form a 90° domain wall at the interface 
between the bands (figure 9). Furthermore, we require full compatibility of 
the spontaneous strain and polarization, which means that the 180° domain 
wall of the two bands meet in one point and the volume fractions of the 180° 
domains within the bands are equal. 

In fact, the question which type of domains are forming the fine structure, and 
which combination forms the bands, is not important in our case, as there is 
no length scale in our model. For intrinsic linear properties of the domain 
structure, the only relevant additional properties are the volume fractions. 
Figure 9 shows two equivalent examples of hierarchical domain structures 
consisting of combined 180° and 90° domains. 

Starting from the 180° structure, described in sect. 4.2, the orientation of one 
band is obtained by a rotation by « 45° about the x<i axis. Then, the properties 
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Fig. 9. Two equivalent examples of hierarchical domain structures 
of the laminar structure from 4.2 have to be rotated using the rotation matrix 
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The properties of band B are obtained applying orientation relationship ac- 
cordingly (22) and with (23). Depending on the volume fractions, the average 
polarization varies within the square spanned by the four individual polariza- 
tion vectors. The average spontaneous strain is reduced in x\ and X3 direction, 
but remains unchanged in perpendicular direction. If a crystal is embedded 
within a elastic neighborhood, hierarchical micr o-structures can reduce the in- 
ternal stresses and the internally stored energy fjArltL 119901 ) . In this case, such 
an in-plane arra ngement of domains d oes not permit an optimal relaxation of 
internal energy ( Arlt and Saskol . Kl 



A domain arrangement, which is often observed in ceramic materials, may be 
characterized by a 3-dimensional arrangement of ferroelectric domains (see fig. 
10). Two bands of 90° domain structures ([001]/[100] with (101) interface and 
[001]/ [010] with (011) interface) are combined in such a way, that the interface 
between the bands is the pseudo-cubic (110) plane. Two domains form a 90° 
domain structure across the band interface ([001]/[00l]) and two domains form 
a 90° domain structure ([100] /[0l0]). Assuming equal volume fraction for all 
four domain, then the first domain band has an average polarization along 
[101] and the second band along [Oil] resulting in an average polarization of 
the hierarchical structure P s * = (1, — l,0)P s /4. 



Starting from the 90° domain structure from sect. 4.1, a series of three rota- 
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Fig. 10. Three dimensional arrangement of four dif ferent polarization dir ections, 
with effective polarization along pseudo-cubic [110], ( Arlt and Saskol . 1980h 



tions is necessary to transform the crystal into a suitable orientation with the 
interface normal of the (110) plane parallel to the macroscopic x% direction: 
a first rotation by 45° about X2 brings the crystal back into an orientation 
with the pseudo-cubic axes parallel to the axes of the macroscopic coordi- 
nate system. A rotation by -45° about 23 and by 90° about X\ gives the final 
orientation of the band A with 
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The properties of band B are obtained applying orientation relationship ac- 
cordingly (22) and with (23). 

A closely related example is that of a engineered domain structure in a tetrag- 
onal ferroelectric, poled along (111) (fig. 11). The three orientations [100], 
[010], and [001] are arranged in two bands forming a mechanically compatible 
domain structure, but with partially incompatible polarizations on the (110) 
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interface between the band. Band A and its orientation in the macroscopic 
coordinate system is identical to that of the preceding example. Band B con- 
sists of [010] and [001] domains with the domain wall (Oil). The orientation of 
band B in the macroscopic coordinate system is then obtained by applying a 
rotation by -90° about Xi- That part of the interface, where the [100]-domains 
of Band and the [010]-domains of band B meet, is charged. If the volume frac- 
tion of the [001] domains in the two bands is 1/3 and the volume fractions of 
the two bands are equal, then the overall polarization is along [111]. 




Fig. 11. Arrangement of three different polarization directions in a tetragonal fer- 
roelectric, which form an engineered domain structure with polarization along 
pseudo-cubic [111] 



A final example demonstrates a hierarchical domain structure for rhombohe- 
dral crystals poled in the pseudo-cubic [001] direction (fig. 12). It consists of 
two domain bands with 109° domain structures, [111]/[111] with (101) do- 
main wall and [111]/[111] with (101) domain walls. The interface between the 
bands is the (001) plane. Within each band, the polarizations lie with the 
(101) and (101) plane, respectively. For analysis of the effective properties of 
this structure, it is suitable to depart from the 109° domain structure in sect. 
4.3. This has to be rotated by 90° about x 3 and by 45° about x 2 in order to 
arrive at the orientation of the lower half of figure 12. The rotation matrix, 
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which transforms the 109° laminate into band A is 
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The upper half, band B, is again obtained by the orientation relationship of 
a fully compatible structure, (22) and (23). Across the band interface, the 
domains form 71° domain structures with polarizations in the (110) and (110) 
plane, respectively. The average polarization of the composite structure is 
along [001], provided that all domains have the same volume fraction. 




Fig. 12. Rhombohedral domain structure with average polarization along 
pseudo-cubic [001] 
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7 Numerical results 



The theory presented in this paper has been used to evaluate the effective prop- 
erties for various arrangements of ferroelectric domains in tetragonal Barium 
Titanate. This material is chosen as it is well characterized and may serve 
as a prototype material for tetragonal ferroelectrics. The results for distinct 
cases are compiled i n table 1. The single crystal properties are taken from 



(jZgonik et all 119941 ) and ere described in an orientation where [001] ||x3. 



The case of 90° laminate refers to an simple laminar domain structure with 
n 1 1 [101] ||x3 and a volume fraction 50% for each domain type. Additionally, 
the bounds on the effective properties, (Q) and (R), are given. The depen- 
dence of the piezoelectric coefficients d 31 , d^ 2 , dg 3 , d* 2A and d\ b on the volume 
fraction is plotted in fig. 13. In this case, the piezoelectric constants d*^ and d% 3 
increase by more than 400% and 100%, respectively, compared to the single 
crystal values and are independent upon the volume fraction. d% 2 is slightly 
smaller than the single crystal value. These effects are, simple averaging ef- 
fects, i.e. the increase of the d^ ± and d^ 3 are only a result of the rotation of the 
crystal and thus caused by the contribution of the large shear coefficient di$. 



More interesting are the piezoelectric coefficients d 24 and d\ 5 , as they show a 
dependence on the interaction between the domains, i.e. these effective proper- 
ties differ from the lower bound (Reuss approximation). They show a quadratic 
dependency upon the volume fraction with a minimum or maximum, respec- 
tively, at £ = 0.5. The effective d\ 5 is slightly boosted by domain interaction, 
but d 2i is strongly restrained. The elastic constants S 1 ^, and S§ 5 , as well as 
the dielectric constants e* x and e 22 are plotted in fig. 14, which shows, that 
the elastic coefficients are much less effected by domain interactions compared 
with the dielectric coefficients. 



The overall symmetry of a simple laminar domain structure is clearly or- 
thorhombic (point group mm2) with special relationships S12 = S23 and 
5*11 = S33 ■ Obviously, the shear compone nts S44, Sm and Sm a re all differ- 
ent, which is in contrast to the finding of ( Erhart and Caol 1999 ). 



Comparing the effective properties with the upper and lower bounds, one has 
to note, that a particular component not necessarily has to be between the 
values of the bounds, but the differences between the property matrices have 
to be positive semi-definite, which is true in the studied cases. 



Furthermore, the effective properties of a 90° domain structure with charged 
domain walls are calculated. The results obtained from the general formula in 
sect. 3 are calculated with a domain wall normal n||^3 and the symmetry 
relationship (35) and a volume fraction £ = 0.5. The average polarization is 
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parallel to the domain wall. For sake of comparison we have finally rotated the 
results so that the average polarization is in x 3 direction, i.e. the interface nor- 
mal is n = (1,0,0)*. This configuration is characterized by much decreased 
piezoelectric constants, except d^ 2 - It is remarkable, that the sign of is 
changed. In table 2 the effective strains and polarization of the two differ- 
ent 90° domain configurations are compared with the single crystal values. 
It demonstrates, how much the effective polarization is increased due to the 
internal electric fields in the case of a charged domain wall. 



Finally, a comparison of the effective properties of three different hierarchical 
domain structures in tetragonal BaTiC>3 is given in table 3. The properties 
are given with respect to the coordinate system, in which they are calcu- 
lated, i.e. the macroscopic £3 axis is perpendicular to the band interfaces. In 
the first case of a 180°/90° domain strcuture with equal volume fractions, the 
macroscopic piezoelectric effect vanishes. Elastic and dielectric properties have 
othorhombic symmetry mm2. In the following cases o f a fully compat i ble 3- 
dimensional arrangement of 90° domains according to Arlt and Saskol ( 1980h 
and the charged 3-dimensional arrangement of 90° domains, the over-all sym- 
metry is not further investigated. It is however obvious, that the symmetry is 
strongly reduced. In the case of the charged arrangement, we have calculated 
a polarization of P s = (0.2,-0.04,0)', which is not the pseudo-cubic [111] 
direction. This effect is again caused by the large internal electric fields, which 
amplify the dielectric displacements in the macroscopic X\ direction. 
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Table 1 

Linear material properties for barium titanate in different domain configurations: 
Single domain single crystals (s. cryst.), uncoupled domains I and II (rotated single 
crystal) , compatible 90° domain structure and a 90° domain structure with charged 
domain wall, both for a volume fraction £ = 0.5. Elastic constants are the stiffness 
at zero electric field in 10~ 3 GPa -1 , the piezoelectric constants are given in pmV -1 
and the dielectric constants are at zero stress, given in 10 -9 Fm 
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Fig. 13. Piezoelectric coefficients for a 90° domain structure of tetragonal BaTiC>3, 
dependent upon the volume fraction of domain I £ 
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Table 2 

Spontaneous strains and polarizations of tetragonal BaTi03 single crystals, in the 
reference orientation and rotated, and effective strains and polarization of fully com- 
patible and charged 90° domain structures. Volume fraction £ = 0.5. Polarizations 
are given in C/m 2 
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Table 3 

Linear material properties for barium titanate in different hierarchical domain 
figurations: Units as in tab. 1 
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8 Summary 



The effective properties of piezoelectric laminates have been analyzed. The 
analysis is based on the calculation of internal fields and makes use of a simple 
matrix manipulation method. The results are expressed in a compact nota- 
tion which is convenient for numerical implementation and at the same time 
suitable for further analytical treatments. 

A more detailed analysis of fully compatible ferroelectric domain structures 
shows, that the results for arbitrary piezoelectric laminates can be further 
simplified and specific property relationships for rank-1 laminates of tetragonal 
and rhombohedral crystals are derived. In particular, it is shown that only 
the components S44, S45 and S55 as well as en, £i 2 and e 22 are effected by 
interactions between the different domains. All other components can be also 
obtained by simple volume averaging of the single domain properties, i.e. they 
are equal to the Reuss lower bound on the effective properties. 

In contrast, in the case of charged domain structure, where the spontaneous 
strains are compatible but not the polarizations, the conditions for a simplified 
treatment are not fulfilled. 

The method is finally applied to the analysis of various hierarchical domain 
structures (rank-2 laminates, or laminates of laminates). Detailed orientation 
relationships between the particular domains in some important domain pat- 
tern are given to make these structures accessible for the presented method. 
Some numerical results for barium titanate are given in order to illustrate the 
effects of different domain arrangements on the effective properties. 
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A Detailed solution of the internal fields for the arbitrary laminate 



A detailed solution of the internal fields (13) with the conditions (17.a-17.c) 
for the arbitrary laminate is de monstrated in t he following. The solution is 
analogous to that one derived by Kreherl (jl998| ). but stress and electric field 
are chosen to be independent variables. 



Although we have given the result for the interaction matrices and vectors in 
eqns. (18.a-18.e) in an index notation, we will derive the general solution in 
this appendix in a symbolic way, which enables a very lucid derivation of the 
results and is convenient for numeri cal implementation . Th is represent a tion o f 
the problem is equivalent to that of Liu and Li ( 2003| ) and Li and Liu ( 2004 ). 
However, their derivation of the results is slightly different. The results can 
easily be converted in the usual notation, simply by rearranging the vector 
and matrix component in an appropriate way. 



We rearrange the components of the vectors p and q , so that the first four 
components represent the field components, for which the strain and dielec- 
tric displacement are continuous across the interface. For the remaining five 
components, the stress and electric field is continuous. 
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Then, the constitutive law takes the form: 
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Q ba = Q 



ab 1 



^ Q-33 Q34 Q35 Q37 Q38 ^ 

Q34 Q44 Q45 Qn Q48 
Q = Q35 Q45 Q55 Q57 Q58 
Q37 Q34 Q57 Q77 Q78 

^Q38 Q48 Q58 Q78 Q88 j 

Correspondingly, the local fields p( fc ) are expressed by 




The conditions (17.a -17. c) are 

a(k) 



a(l) 



P &A(k) =o 



(A.2) 



= q 

p bik) = p b 
P a(fc) ) = p a 



(17.b) or (A.3.b) yields directly: 

pbb( k ) _ jbb j^baik) _ q 

where I bb = 5y, i,j = 1...5. This results are equivalent to (18. a) and 
(A. 3. a) can be written in the following way, using (A.l) 

a(k) _ Qaa(k)^a(k) _|_ Qab( k ) ^b( k ) _|_ qSa(fc) _ o(l) 

With the abbreviation = q sa(fc) + Q^( k ) p b( k ) we get 



(A.3.a) 
(A.3.b) 
(A.3.c) 



I8.d). 



QOO Wp 



(k) a(k) 



q 



a(l) 



,(*) 



(A.4) 



(A.5) 



p a(fc) = (q-M)^ ( q a(D _ a W) (A 6) 

and with (A.3.c) and the definition 1 R aa(t:) = ^Q aa(7c ^ this gives: 



R aa (q a 



(i) _ n (fc) 



R aa \q a(1) - (R aa a W 



q a(l) = / R aa\ 1 + / R aa\ 1 /r^W 



1 Note, that R aa(fc) represents local properties. The suffix (k) will be dropped in 
the following. 
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We define Q aa = ^R aa ^ and write 

q a(l) = Q<*a p a + gaa /jjaa^k) 

which we insert into (A. 6): 

p a(fc) = paa {Qaa-a + gaa /]jaa a (fc)\ _ & (fc) 



After rearranging this, we find 



(A.7) 



p a(fc) = K««Q^ p » + R<^Q<™ (R aa Q afc(fe) \ p 6 - R aa Q ab{k) p b 



R aa C^ aa /R aa q sa ^\ R aa q sa 



(A.8) 



By comparison with (A. 2) we find the following expressions for A aa , A ab and 



j^aa R, aa Q aa 

n A 

P 



R aa Q 



ab( k ) 



(A.9.a) 
(A.9.b) 
(A.9.c) 



Equation (A. 9. a) is equivalent to (18. b), equation (A.9.b) is equivalent to 
(18. c) and equation (A.9.c) is equivalent to (18. e). 

Using this notation we can also express the effective properties in a symbolic 
way: 



Qaa* Qab* \ 
Qba* Qbb* J 
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K Q ba(k) Q bb(k) 
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= ^Q 6a A afe - 


hQ bb ) 




= (Q ba R aa \ 
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And the effective spontaneous strain and polarization are: 

q ~ V q s * b J ~ H Q6a(fc) Q66(fc) / V pAKfe) / ' + U qS / ' 

q s*a = /Qaa p A»\ + ( q sa^ = gaa ^«a q sa^ 

It seems not to be useful to give explicit equations for the various components 
of this general solution, as the expressions are quite lengthy. We have not as- 
sumed any symmetry relation for the crystal and its property tensors and will 
therefore not expect to find any particular symmetry changes in this results. It 
should be noted that the matrices Q aa ( k \ Q ab ^ and R aa are constants for the 
various phases k. On the contrary, Q act is a over-all property which depends on 
the volume fractions Finally, the interaction matrices A aa ^ and A a& ^ 
hold for a particular phase k, but depend on the volume fraction £( p ) of all 
phases and are not particular properties of the phase k alone. 

We can proof the general solution by showing that the micro-macro relations 
are fulfilled for all components: 

p = (p) = (A) p + (p A ) 

which is true, if (A) = I and (p^ = 0. This can be shown easily. Additionally, 
one can show, that q a are homogeneous, i.e. they depend on macroscopic 
quantities only: 

Homogeneous material properties 

If Q is homogeneous, i.e. all lamellae have the same linear properties, then 
follows: 

gaa = gaa 

A ab = 0, 

where I aa is a 4 x 4 identity matrix, which shows, as expected, that A = I, 
with I beeing a 9 x 9 identity matrix. Then 

Q* = (QA) = Q. 
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If Q is homogeneous, but the spontaneous strain and polarization q s are not 
homogenoeus, then 

paA _ paa r^q so ( fe A — q sa ( fc )^ 

which says, that the local stress and electric field are proportional to the dif- 
ference between the average spontaneous strain and polarization and the local 
one. We find that for homogenous linear properties the effective spontaneous 
strain and polarization are equal to their simple volume average: 

q s * = (q s ) (A. 10) 

Partial homogeneity 

It is interesting, and for the case of ferroelectric domain structures impor- 
tant to note, that considerable simplifications of the general equations for 
the effective properties may be found, even if not all components of Q are 
homogeneous. If only the Q aa , and therefore R aa , are homogeneous, then 



Q aa R aa = I aa and we get following relationships: 

A aa = R aa Q aa =I (A. 11) 

A ab = paa {/q^A _ Q«t( fc )^ (A. 12) 

p aA = R aa ((q sa(fc) ) -q sa(fe) ) . (A.13) 
and for the effective properties 

Q*aa = ^aa^ = Qaa (A. 14) 

Q* ab = (Q ab ) (A.15) 

Q* ba = (Q ba ) = (Q* ab y (A. 16) 

Q* bb = R aa (Q ab ) - (Q ba R aa Q ab ) + (Q bb ) (A. 17) 

q s*a = ^ ^ A18 ^ 



If additionally the components q sa are homogeneous, then the latter equation 
simplies to 




For further discussions it is useful to split the effective linear properties into 
the simple average (Q) and an additional contribution A. It is then convenient 
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to express the interaction matrix by A = I + A': 



Then 



A' = 1 + 



A ab 




Q* = (Q A) = (Q (I + A')) = (Q) + (QA'> 





A = (QA'> = Q* - (Q) 

A bb = /q6o\ paa /qo6\ _ /Q&apaaQ 



A bb 
ba t> aa(~\ab\ 
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